UK/99-20 

Quark Orbital Angular Momentum from Lattice QCD 

N. Mathur a - 6 , S. J. Dong a , K. F. Liu Q ' c , L. Mankiewicz d , and N. C. Mukhopadhyay 6 
a Dept. of Physics and Astronomy, University of Kentucky, Lexington, KY 40506 
b Dept. of Physics, Applied Physics and Astronomy, RPI, Troy, NY 12180 
C SLAC, P. O. Box 4349, Stanford, CA 94309 
d Physics Dept., Technical Univ., Munich, D-85747, Garching, Germany 



ON 

On 

o 

Q 

o 



> : 

On . 
OO ■ 
(N ■ 

<n : 

On ■ 

5 : 

^ : 

Oh! 
I 

^ ■ 



X 



We calculate the quark orbital angular momentum of the nucleon from the quark energy-momentum 
tensor form factors on the lattice. The disconnected insertion is estimated stochastically which 
employs the Z2 noise with an unbiased subtraction. This reduced the error by a factor of 4 with 
negligible overhead. The total quark contribution to the proton spin is found to be 0.30±0.07. From 
this and the quark spin content we deduce the quark orbital angular momentum to be 0.17 ± 0.06 
which is ~ 34% of the proton spin. We further predict that the gluon angular momentum to be 
0.20 ± 0.07, i. e. ~ 40% of the proton spin is due to the glue. 



PACS numbers: 12.38.Gc, 13.88.+e, 14.20.Dh 

The spin content of the proton remains a challeng- 
ing problem in QCD both experimentally and theoret- 
ically [jjj. The surprisingly small contribution from the 
quark spin revealed by the polarized deep inelastic scat- 
tering experiments || (world average: E = 0.25 ± 0.10) 
has stimulated a great deal of interest in the understand- 
ing of this 'proton spin problem'. While the lattice QCD 
calculations |^] confirmed the small quark spin content in 
agreement with experiments, there is little consensus on 
where the rest of the proton spin resides. There have been 
suggestions based on the Bjorken sum rule , the parton 
evolution (5) , the chiral quark model || and skyrmion J?]] 
that the quark orbital angular momentum in the nu- 
cleon can be substantial. It is further proposed that the 
off-forward parton distributions from the deeply virtual 
Compton scattering can be used to measure the quark 
orbital angular momentum distribution and thereby its 
moments B. 

In this letter, we shall report the first lattice calculation 
of the quark energy-momentum tensor form factors which 
admits the extraction of the total quark angular momen- 
tum. Combining the lattice calculation of the quark spin 
content H , we obtain the quark orbital angular momen- 
tum and thereby predict the gluon angular momentum in 
the nucleon from the spin sum rule. It turns out that the 
quark orbital angular momentum is indeed quite large. It 
constitutes ~ 35 % of the proton spin and the gluon an- 
gular momentum is predicted to make up the remaining 
~ 40% of the proton spin. 

It has been shown recently || that the total angular 
momentum in QCD can be decomposed into three pieces 
in a gauge invariant way 



J= I d 3 x ^ V'TTsV' + / d 3 x ^ji{x x (—iD)}i/) 



d 3 x [x x (E x B)] = -E + L q + J g , (1) 



The quark spin contribution ^E is defined through the 
forward matrix element of the flavor-singlet axial-vector 
current. Similarly, L q , the quark orbital angular momen- 
tum, is defined from the operator L q and J g , the gluon 
angular momentum, is defined from J g which is asso- 
caited with the gluon Poynting vector E x B. The total 
quark angular momentum is J q = iE + L q . Whereas, 
the total gluon angular momentum J g cannot be further 
decomposed into gluon spin and gluon orbital angular 
momentum without explicit gauge dependence. 

J q can be obtained from the energy-momentum tensor 
form factors || . One can write the gauge invariant quark- 
gluon energy-momentum tensor as 
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= \^[D - "^D] v )ip + F„ a F va - ^S^ V F 2 , (2) 

where the first part is the quark energy-momentum ten- 
sor and the second one is that of the gluon. Form factors 
of the energy-momentum tensor for either the quark or 
the gluon are defined by M, 



< P, s\T^(0)\p', s' >= u(p, s)[Ti(g 2 ) i^$> v 



T 2 (q 2 )p( tJ .i<7u) a qa/2'rn N H T 3 (q 2 )(q t _ l q u - S^q 2 ) 



m N T±(q 2 )5 fl v]u(p ', s'), 



(3) 



where p^ = (p^+p' ft )/2, q^^p^-p'^ and u(p) is the nu- 
cleon spinor. It can be proven that, for polarized target, 
the total angular momentum of the quarks or gluons 
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Therefore, the total angular momentum of the 
quarks/gluons can be calculated at the q 2 — > limit of 
the form factor l/2[T?' 9 {q 2 ) + T%' 9 (<?)}. 

In Eq. (H), the energy momentum tensor is expressed 
in terms of four form factors; whereas, for calculating the 
angular momentum we need only two, namely T\ and T 2 . 
Tij, with j in the 3-direction, does not admit the T4 term. 
To remove T3, we choose the momentum transfer to be 
orthogonal to the j direction. In order to get the re- 
quired Ti(q 2 )+T 2 (q 2 ), we calculate the three point func- 
tion Gmt^ -.n (*2 1 t\i Pi —Q) for the operator T^. The three 
point function has two parts: connected insertion (CI), 
due to the valence and cloud quarks, and disconnected 
insertion (DI) attributed to the sea quarks M (Fig.l). 
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FIG. 1. Quark skeleton diagrams for the three-point func- 
tion GNT 4j N{t2,ti,p, —q). (a) is the connected insertion, (b) 
is the disconnected insertion. 



For CI, we express the observable T(q 2 ) = ± [Ti(g 2 ) + 
T2(q 2 )] in terms of the ratio of three- and two-point func- 
tions : 
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where T m and T e are the spin polarized and unpolar- 
ized projection operators pfj[ . From the above ratio, we 
calculate the lattice Tci(q I )aX different q 2 and then ex- 
trapolate them to q 2 — » limit to obtain the CI part 
of the lattice quark total angular momentum Jq^ci- Re- 
sults are obtained for relatively light Wilson quarks with 
k = 0.148,0.152,0.154 and 0.155 which correspond to 
quark masses of about 376,210, 124 and 80 MeV re- 
spectively. Fig. 2 shows the dipole fitting of Tci(q 2 ) at 
different q 2 . Following the calculation for the point- 
split Wilson current JlOj ] , the tadpole improvement fac- 
tor l/8K c (±TrU p i aq )^(K c = 0.15684) is included in 
the unrenormalized J q ,ci- The chiral limit for J q ,ci 
is taken from a linear dependence on the quark mass 
m q a for these four k values (see Fig. 3). To account for 
the correlations, both the dipole fitting and the chiral 
extrapolation are done with the covariance matrix and 
the final error at the chiral limit is obtained from the 
jackknife procedure. The dipole mass is found to be 
0.88 ± 0.07 GeV. Finally, to obtain the result in the ~MS 



scheme at 1/a = 1.74GeV, we multiply the J q ,ci by the 
tadpole improved renormalization constant for the oper- 
ator Tij which has been calculated [|lj perturbatively to 
be Z = 1.045 and obtain the CI part of the quark an- 
gular momentum J q ^ C i = 0.44 ± 0.07. Being ~ 90% of 
the nucleon spin, this almost saturates the spin sum rule. 
We also performed monopole fitting for the Tci(q 2 ) and 
found an order of magnitude larger \ 2 ■ The calculation 
is done on a quenched 16 3 x 24 lattice at (3 — 6.0 with 
Wilson fermions for 100 configurations. 
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FIG. 2. Dipole fitting for Tci{q 2 ) at 4 different k values. 
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FIG. 3. Chiral extrapolation for J q ,ci as a function of the 
quark mass. The value at the chiral limit is indicated by •. 

For the DI contribution, we follow the calculations for 
the flavor-singlet axial coupling ||, the irNa term [ EjL 
and the strangeness magnetic moment jl3j by summing 
the current insertion time t\ from the nucleon source to 
its sink in the corresponding ratio in Eq. (^|) for the DI 
(Fig. 1(b)) to gain statistics. In this case, the ratio leads 
to const + l/2€jk m qkTDi(q 2 ) h- We take the average of 
the 3 polarization directions. The DI result is calculated 
from the slope of this summed ratio with respect to t2- 
To evaluate the trace of the quark loop in Fig. 1(b), we 
adopt the same stochastic algorithm with the Z2 noise 
estimator |l4j as in other DI calculations P, pll| , |l3| ] . In 
addition, we shall use two more techniques to reduce the 
errors from the stochastic algorithm. First one is to ob- 
serve that from the charge conjugation and Euclidean 
hermiticity (CH symmetry), the three-point function for 
the DI which is the product of the quark loop and the 
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nucleon propagator is real. On the other hand, the Eu- 
clidean hermiticy itself dictates that the loop should also 
be real. Therefore, we need only multiply the real part 
of the loop with the real part of the nucleon propagator 
and neglect the product of the imaginary parts of them 
which only introduces noise to the signal. Secondly, we 
employ an unbiased subtraction method which has been 
used in the calculation of the fermion determinant [fl5[ . 
The trace of the inverse matrix A^ 1 can be estimated 
stochastically as follows 



Tr(A- v ) = E[< vHA- 1 - A * >1> 



(6) 



where 77's are Z2 noise vectors, are a set of P 

traceless matrices and A,'s are the variational parame- 
ters which are determined by reducing the variance of 
the three-point function over the gauge configurations. 
In practice, we found that a judicious choice of 0^> is a 
set of traceless matrices from the hoping expansion of the 
propagator. Since they match the off-diagonal behavior 
of the matrix A -1 , they can offset the off-diagonal con- 
tribution to the variance p5|. This method proves to 
be very efficient in reducing the error of the DI calcu- 
lation with negligible overhead. After implementing the 
CH and H symmetries and the unbiased subtraction with 
traceless matrices obtained from just the first two terms 
of the hopping expansion, we obtain a reduction of error 
of 3 - 4 times. Fig. 4 shows a plot of the summed three- 
point to two-point function ratio (as in Eq. (j^ for the 
DI) vs. ti for k = 0.154 and \q\ — 2ir/La with and with- 
out subtraction. One can see that the error bars before 
subtraction are much larger and only with the subtrac- 
tion does one get a reasonably good slope as illustrated by 
the fitted straight line. With the help of this subtraction 
procedure we calculate slopes for other q 2 and for other 
k's which are shown in Fig. 5. We use fixed source and 
vary the sink position t%. From ti = 8 on, the nucleon 
becomes isolated from its excited states j[(J. Hence, the 
slopes are fitted in the region ti > 8 (Fig. 5). Next, 
Tr>i{q 2 ) is fitted with a monopole form in q 2 as in the 
other DI calculations |l2|]l5| ]. These are plotted in Fig. 
6. Similar to the CI case, we also use covariant matrix 
fitting and the final error bars are obtained by the jack- 
knife method. A finite mass correction factor from the 
triangle diagram |lrj| ] is introduced while extrapolating 
to the chiral limit with a linear m q a dependence. This 
is shown in Fig. 7(a). The strange quark contribution 
is obtained by fixing the sea quark mass at k s = 0.154 
and extrapolate the valence k v from 0.148,0.152,0.154 
to k c . This is shown in Fig. 7(b). It is interesting to 
point out that the result is fairly independent of the sea 
quark mass in the fermion loop ( Fig. 1(b)). Comparing 
Fig. 7(a) where the valence and sea quarks are kept the 
same (k v = k s ) and Fig. 7(b), we see that albeit the sea 
quark mass in Fig. 7(a) changes by a factor of 3, the 
result still coincide with those in Fig. 7(b) for each of 
the valence-quark case. This shows that the DI depends 



on the valence-quark mass but is almost independent of 
the sea-quark mass. 
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FIG. 4. The summed ratios of Eq.(|^) for the DI are plotted 
for different time slice t2 with and without unbiased subtrac- 
tion. Ratios without subtraction are shifted slightly towards 
right. 
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FIG. 5. The summed ratios of Eq.(fe|) for the DI are plotted 
for |q| = 2n/La. 
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FIG. 6. Monopole fitting for T D i(q 2 ). Tb/(0) values are 
obtained by extrapolating to q 2 — > 0. 
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FIG. 7. (a) The ^faitice J q ,Di as a functiJn'W the quark 
mass m q a. The quark masses in the valence and the sea are 
kept the same, (b) J a vs the valence quark mass with the 
sea quark mass fixed at k s = 0.154. The chiral limit value 
(k v = At c ) is indicated by •. 
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This is reminiscent of the sea-quark contribution in 
the flavor-singlet g° A calculation S where the sea-flavor 
independence was first observed. 

The breakdown of the quark angular momentum J q 
into the quark spin and the quark orbital angular 
momentum is given in Table f . From the CI calculation 
we obtain the valence and cloud quark contributions to 
the quark angular momentum J q ,ci = 0.44 ± 0.07 which 
is ~ 90% of the total proton spin and it almost saturates 
the spin sum rule in Eq. (|l|) by itself. A previous cal- 
culation of the flavor-singlet axial current on the same 
set of lattices shows that \T>ci — 0.31 ± 0.04 Since 
Jq,ci — ^ci + Lq t ci, we obtain the CI part of the quark 
orbital angular momentum L Qt ci — 0.13 ± 0.07. Thus, 
for valence and cloud quarks, about 70% of J q .ci comes 
from the quark spin and the 30% is due to the orbital an- 
gular momentum. From the DI calculation, we find that 
the total quark angular momentum J q ,Di, like the quark 
spin ^Tijji, is also flavor symmetric within errors. In 
fact, J u ,di, Jd,Di, and J s are all equal to —0.047 ±0.012. 
Together, the total DI is J q ^ DI = -0.14 ±0.04. Subtract- 
ing the DI of the quark spin = —0.18 ± 0.03 from 
Jq,Di, we obtain the orbital angular momentum contri- 
bution from the sea quarks to be L q ^i = 0.041 ± 0.035. 
It is interesting to note that it is consistent with zero 
with a central value which is a factor of 4.5 smaller 
than the spin content of the sea quarks. Adding CI and 
DI contributions together, we obtain J q = 0.30 ± 0.07 
and, thus, we predict the gluon angular momentum 
J g = 0.5 - 0.30 ± 0.07 = 0.20 ± 0.07 from the spin sum 
rule (Eq. (0)). 

To conclude, the total angular momentum of the 
quarks is calculated to be J q = 0.30 ± 0.07, i. e. ~ 60% 
of the proton spin is attributable to the quarks. Since 
the quark spin content is calculated previously to be 
iE = 0.13±0.06 ||, we obtain the quark orbital angular 
momentum L q = 0.17 ± 0.06. Therefore, about 25% of 
the proton spin originates from the quark spin and about 
35% comes from the quark orbital angular momentum. 
The gluon angular momentum contribution is predicted 
from the spin sum rule to be J g — 0.20 ± 0.07, i. e. ~ 
40% of the proton spin is due to the glue. Since the or- 
bital angular momentum of the sea quarks turns out to 
be quite small, the sea flavor independence of J q ,Di re- 
confirms the sea flavor independence of the quark spin, 
namely Au(DI) = Ad(DI) ~ As observed in the previ- 
ous lattice calculation j|. 



TABLE I. The Breakup of Quark Angular Momentum 





J q 


?^ 


L q 


u + d{CI) 


0.44(7) 


0.31(4) 


0.13(7) 


u/d(DI) 


-0.047(12) 


-0.062(6) 


0.015(12) 


s 


-0.047(12) 


-0.058(6) 


0.011(12) 


u + d + s(DI) 


-0.14(4) 


-0.18(3) 


0.041(36) 


Total 


0.30(7) 


0.13(6) 


0.17(6) 



In addition, the fact that J q .ci almost saturates the spin 
sum rule and the sea quark orbital angular momentum is 
small, the gluon angular momentum and the sea quark 
spin largely cancel each other. Both the flavor indepen- 
dence in the sea and the cancellation between the gluon 
angular momentum and the sea quark spin suggest that 
it is related to the anomaly and anomalous chiral Ward 
identity. 

It is important for the future experiments to measure 
the quark orbital angular momentum and the gluon angu- 
lar momentum in order to conclude the study of the spin 
structure and content of the proton. The present work 
is based on the quenched approximation and is subject 
to the large volume, finite lattice spacing and chiral ex- 
trapolation corrections. We shall address these issues in 
a future study with the overlap fermion |l7| which has 
nice chiral and scaling properties. 
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